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ABSTRACT 

Using the Tibet-Ill air shower array, we search for TeV 7-rays from 27 poten- 
tial Galactic sources in the early list of bright sources obtained by the Fermi Large 
Area Telescope at energies above 100 MeV. Among them, we observe 7 sources 
instead of the expected 0.61 sources at a significance of 2a or more excess. The 
chance probability from Poisson statistics would be estimated to be 3.8xl0~ 6 . 
If the excess distribution observed by the Tibet-Ill array has a density gradient 
toward the Galactic plane, the expected number of sources may be enhanced 
in chance association. Then, the chance probability rises slightly, to 1.2xl0~ 5 , 
based on a simple Monte Carlo simulation. These low chance probabilities clearly 
show that the Fermi bright Galactic sources have statistically significant correla- 
tions with TeV 7-ray excesses. We also find that all 7 sources are associated with 
pulsars, and 6 of them are coincident with sources detected by the Milagro exper- 
iment at a significance of 3a or more at the representative energy of 35 TeV. The 
significance maps observed by the Tibet-Ill air shower array around the Fermi 
sources, which are coincident with the Milagro >3a sources, are consistent with 
the Milagro observations. This is the first result of the northern sky survey of 
the Fermi bright Galactic sources in the TeV region. 

Subject headings: gamma rays: observations — pulsars: general — supernova 
remnants 
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1. Introduction 

The Fermi Gamma-ray Space Telescope (Fermi), succeeding the Energetic Gamma Ray 
Experiment (EGRET), was launched in June 2008 to cover the energy range of 20 MeV 
to 300 GeV, with a sensitivity approximately a hundred times better than that of the 
EGRET. The Large Area Telescope (LAT) on board the Fermi surveyed the entire sky for 
3 months, after which the 205 most significant sources were published in a bright source 
list above 100 MeV at a significance greater than ~10er (lAbdo et al. 2009al ). Remarkably, 
this survey detected many new 7-ray pulsars. A typical 95% uncertainty radius of source 
position in this list is approximately 10' and the maximum is 20'; these values are greatly 
improved compared to those of the EGRET. This provides a more accurate, unbiased 
search for common sources across multi wavelengths, compared with the EGRET era. 
Recently, the Milagro experiment observed 14 of the 34 Fermi sources selected from the 
l ist at a false-positi ve significance of 3a or more at the representative energy of 35 TeV 



f lAbdo et al. 2009bh 



In this paper, we report on a search for TeV 7-ray sources in the Fermi bright source 
list with the Tibet-Ill air shower array (Tibet-Ill array). We also discuss simple statistical 
tests for our results and possible coincidences with the Milagro observations. 



2. Tibet-Ill Air Shower Array 

The Tibet air shower array has been operating at Yangbajing Cosmic Ray Observatory 
(90?522 east, 30?102 no rth; 4300 m above sea level) in Tibet, China, since 1990 



( lAmenomori et al. 19921 ). We observe cosmic rays and 7-rays using the extensive air 
shower technique of a scintillation detector array with a duty cycle of about 24 hours every 
day regardless of weather conditions, and with a wide field of view of about 2 sr. These 
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capabilities give us an unbiased survey of the northern sky. After several upgrades, the 
Tibet-Ill array was completed and started collecting data in late 1999. This array consists 
of 533 plastic scintillation detectors of 0.5 m 2 placed at grid poin t 7.5 m apart, and its 
coverage area is approximately 22,050 m 2 (jAmenomori et al. 20031 ). Each detector, called a 
fast-timing (FT) detector, has an FT photomultiplier tube to collect scintillation photons. 
The number of air shower particles and the arrival timing of particles at each detector 
are recorded, allowing us to estimate primary cosmic ray or 7-ray direction and energy 
for each air shower. The systematic uncertainty of the absolute energy scale observed by 
the Tibet-Ill array in the m ulti-TeV region is calibr ated to be less than ±12% using the 



Moon's shadow observation (jAmenomori et al. 20091 ) . The single event angular resolution 



is estimated to be 0?9 for modal energy 3 TeV, although it depends on the observed number 
of air shower particles. The systematic pointing error is also estimated to be smaller than 



0?011. These are verified by the Moon's shadow observation (jAmenomori et al. 20091 ). 



Using the Tibet-Ill and previo us arrays, we have succ e ssfully observed TeV 7- ray 



sources, s uch as the Crab Nebul a (jAmenomori et al. 1999 



Mrk 501 (jAmenomori et al. 20001 ). and Mrk 421 (jAmenomori et al. 20031 ). We have also 



Amenomori et al. 20091 ) 



successfully drawn a 



precise two-dimensiona 



map of the large-scale cosmic-ray anisotropy 



in the northern sky (jAmenomori et al. 20061 ). where we first pointed out new small-area 



enhancements in the Cygnus arm direction at multi-TeV energies. One of the enhancements 
is coincident with MGRO J2019+ 37, which was est ablished recently by the Milagro 



experiment as a TeV 7-ray source (jAbdo et al. 20071). It is worth noting that the Tibet 



AS7 experiment has reported several times (jZhang 2003 



Zhang et al. 2005 



Zhang 20051 ) 



on the marginal excesses from the direction closing to MGRO J1908+06/HESS J1908+063 
before the final discovery made by the Milagro experiment. 
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3. Air Shower Data Analysis 



We analyze the air shower dataset collected by the Tibet-Ill array during 1915.5 live 
days from 1999 November through 2008 December. To extract an excess of TeV 7-ray 
air shower events coming from the direction of a target source in this analysis, we adopt 
almost the sam e event selections and the background estimatio n method published in our 



previous work (jAmenomori et al. 2003 



Amenomori et al. 20091 ). We use air shower events 



with ^2 Pft > 10 1 ' 25 as the primary energy reference, where the size Pft is defined 
as the sum of the number of particles per m 2 for each FT detector. The modal 7-ray 
energy, assuming the Crab's orbit and integral spectral index —1.6, is calculated to be 
approximately 3 TeV by the Monte Carlo simulation. The modal 7-ray energy depends on 
the declination and is estimated to be ~3 TeV for a declination band from 20° to 40° and 
~6 TeV for declinations at 0° and 60°, respectively. The search window radius centered 
at the target source is expressed by R(%2 Pft) = 6.9/ a/S Pft degrees, which is shown 
to maximize the signal-to-noise (S/N ) ratio by Monte Carlo study assuming a point-like 
7-ray source (jAmenomori et al. 2003 ). Therefore, an excess might be underestimated if the 
target source actually extends beyond our angular resolution size. 

The target sources in the Fermi bright source list are chosen as confirmed or 
potential Galactic s ources in a similar way to that employed by the Milagro observation 



( lAbdo et al. 2009bl ). Out of the 205 most signific ant sources in the Fermi bright source 
list, 83 are not identified as extragalactic sources (lAbdo et al. 2009al ). Among these 83, we 
select 27 sources in the declination band between 0° and 60°, corresponding to our sensitive 



field of view for TeV 7-ray sources. 
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4. Results and Discussion 

The Tibet-Ill array observation of the selected 27 Fermi bright Galactic sources is 
summarized in Table [U where 13 of the selected sources are classified as pulsars (PSR), 5 are 
supernova remnants (SNR), and 9 remain unidentified but are pote ntial Galactic sour ces, 



as they are mostly concentrated in the Galactic plane <~ 20°) (lAbdo et al. 2009al ). As 
a result of this excess search for these sources, we find no statistically significant evidence 
for TeV 7-rays from other individual sources except for the Crab, which is recognized as the 
brightest standard TeV source. 

Subsequently, the distribution of the observed significance is examined for statistical 
consistency with the normal Gaussian. Figure [T] shows the significance distribution of the 
27 sources observed by the Tibet-Ill array. One can see that the 7-rays from the Crab are 
detected at a sufficiently high significance of 6.9a. It should be emphasized that we observe 
7 sources including the Crab at a significance of 2a or more in this distribution, against an 
expected 0.61 sources (upper probability of 2a multiplied by 27 sources) from the normal 
Gaussian. The chance probability from Poisson statistics would be estimated as 3.8xl0~ 6 . 
With the Crab excluded, the chance probability would be estimated as 3.6xl0 -5 . This low 
chance probability clearly shows that the Fermi bright Galactic sources have statistically 
significant correlations with the TeV 7-ray excesses. In order to check the possible bias of 
the data sample, spatially independent dummy sources are selected from the northern sky 
except for the Crab, Mrk 421, and two famous l arge-scale cosmic-ray an isotropy regions 



known as the Loss-Cone and the Tail-In regions (lAmenomori et al. 20061 ). The Loss-Cone 
and the Tail-In regions are separated from the selected 27 Fermi sources. As a result, the 
significance distribution of ~2000 dummy sources is consistent with the normal Gaussian 
with a mean value of m = —0.010 ± 0.025 and a standard deviation of a = 1.027 ± 0.019. If 
the significance distribution observed by the Tibet-Ill array has a density gradient toward 
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the Galactic plane, the expected number of sources at 2a or more may be enhanced in 
chance association. To check this, we perform a simple Monte Carlo simulation in a similar 
way to that implemented by Romero et al. (1999). We generate 2000 dummy-source lists of 
the 27 Fermi bright Galactic sources, where the Galactic latitude distributions retain the 
form of the actual histograms of the Fermi sources with the Galactic longitudes randomly 
set to new distributions within our field of view. We count the number of >2a sources in 
the Tibet-Ill data according to each dummy-source list, and calculate that the average 
expected number of sources at 2a or more is 0.73±0.02 in the chance association. In this 
case, the chance probability associated with 7 sources or more goes up slightly, to 1.2xl0~ 5 , 
while it becomes 9.7xl0~ 5 with the Crab excluded. 

The Milagro observation found 14 out of 34 Fermi sources at a significance of 3a 
or more, and its sensitivity is approximately two or three times better than that of the 
Tibet-Ill array. Hence, our threshold significance 2a, which corresponds to ~30% of the 
Crab flux assuming a point-like source, should be a reasonable value. We note that the flux 
of OFGL J2020. 8+3649 seems to be quite low compared with the Milagro's, since the flux 
measurement of OFGL J2020.8+3649 with the Milagro experiment is (67±7)% of the Crab 
flux above 35 TeV, while our flux is (30±14)% of the Crab flux above 3 TeV. The statistical 
difference between them is calculated to be 2.3a. This difference may be explained by either 
statistical fluctuations, a harder energy spectrum than that of the Crab, or an extended 
source instead of the assumed point-like source in this analysis. 

We also find that all 7 sources at 2a or more are associated with pulsars, and 6 of 
them are coincident with sources detected by the Milagro experiment at a significance 
of 3a or more at the representative energy of 35 TeV. The remaining source still has 
a positive significance of 1.4a measured by the Milagro experiment. Furthermore, the 
latest Fermi LAT observations detected 16 7-ray pulsars in the blind frequency searches 
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( lAbdo et al. 2009a ). among which only LAT PSR J2238+59 is a new pulsar not included in 
the Fermi bright source list. We also find a 2.5a excess associated with LAT PSR J2238+59, 
as listed in Table [TJ The location of this pu lsar is observed at a significance of 4.7a by 



the Milagro experiment (jAbdo et al. 2009b! ) . In this connection, the first Fermi LAT 



catalog including 4 6 7- ray pulsars has been published using the first 6 months of data 
(jAbdo et al. 2009dl ) r 8 l . We remark tha t no pulsed emis sion has ever been detected from 
7-ray pulsars above a few tens of GeV (lAliu et al. 20081 ). 



Around the Fermi Galactic bright sources with >2a significance by the Tibet-Ill data 
and >3cr by the Milagro data except for the Crab, we compare significance maps between 
the Tibet-Ill array (a)-(d) and the Milagro experiment (a')-(d') taken from Abdo et al. 
(2009b) in Figure [2j Each image has a 5° x 5° region including one or two Fermi sources. 
It is remarkable that the Tibet-Ill array obtains images consistent with those observed in 
the Milagro experiment. Besides, the maximum significance positions obtained by both 
the Tibet-Ill array and the Milagro experiment might be shifted from the pulsar positions. 
In fact, recent imaging air Cherenkov telescopes also discovered many candidates for TeV 
pulsar wind nebulae (PWNe), w hich are displaced within a few tenths of degr ee from 



the pulsars in the southern sky (lAharonian et al. 2006 



Aharonian et al. 20071 1 . Thus 



these observations would imply that the excesses are possible candidates for TeV PWNe. 
The correlation between TeV and GeV 7-rays is being realized by the wide sky survey 



28 Note added in press: in this catalog, 18 pulsars are located in the Tibet-Ill field of 
view. Among them, 14 pulsars are listed in Table [1] including LAT PSR J2238+59. We find 
no significant excess at 2a or more from the other 4 pulsar locations not listed in Table [TJ 
Using the first Fermi LAT pulsar catalog, we find that 8 out of 18 pulsar locations have a 
significance of 2a or more. In this case, the chance probability falls to 1.4xl0 -8 , while it 
becomes 1.8xl0~ 7 with the Crab excluded. 
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instruments, such as the Tibet-Ill array and the Milagro experiment, in the early Fermi era. 

The collaborative experiment of the Tibet Air Shower Arrays has been performed 
under the auspices of the Ministry of Science and Technology of China and the Ministry 
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Science and Technology, by Grants-in-Aid for Science Research from the Japan Society 
for the Promotion of Science (JSPS) in Japan, and by the Grants from the National 
Natural Science Foundation of China, the Chinese Academy of Sciences, and the Ministry 
of Education of China. K. Kawata is supported by Grant-in-Aid for JSPS Fellows 21-9437. 
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Fig. 1. — Histograms show significance distribution of the Fermi bright sources observed by 
the Tibet-Ill array. The dashed curve is the expected normal Gaussian distribution. 
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Fig. 2. — Comparisons of significance maps around the Fermi sources between the Tibet- 
Ill array (a)-(d) and the Milagro experiment (a')-(d') taken from Abdo et al. (2009c). 
Selected are Fermi sources with >2a significance by the Tibet-Ill array and >3cr by the 
Milagro experiment except for the Crab. White points in each image show the Fermi source 
positions: (a) (a 1 ) J1907.5+0602/J1900.0+0356; (b)(b') J0634.0+1745 (Geminga); (c)(c') 
J2021.5+4026/J2032.2+4122; (d)(d') J2020.8+3649. The horizontal axis, vertical axis, and 
color contours indicate the right ascension, declination, and significance, respectively. 
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Table 1. Summary of the Tibet-Ill array observations of the Fermi sources. 



Fermi LAT 








Tibet-Ill 


Milagro a 


Source 


Source 


Class 


R.A. 


Dec. 


Signi. 


Signi. 


Associations 


(OFGL) 




(deg) 


(deg) 


(*) 


(*) 




J0030.3+0450 


PSR 


7.600 


4.848 


1.7 


-1.7 




J0357.5+3205 


PSR b 


59.388 


32.084 


-1.7 


-0.1 




J0534.6+2201 


PSR 


83.653 


22.022 


6.9 


17.2 


Crab 


J0617.4+2234 


SNR 


94.356 


22.568 


0.2 


3.0 


IC 443 


J0631.8+1034 


PSR 


97.955 


10.570 


0.3 


3.7 




J0633.5+0634 


PSR b 


98.387 


6.578 


2.4 


1.4 




J0634.0+1745 


PSR 


98.503 


17.760 


2.2 


3.5 


Geminga 


J0643.2+0858 




100.823 


8.983 


-1.2 


0.3 




J1830.3+0617 




277.583 


6.287 


-0.2 


0.2 




J1836.2+5924 


PSR b 


279.056 


59.406 


-0.3 


-0.9 




J1855.9+0126 


SNR 


283.985 


1.435 


0.7 


2.2 


W44 


J1900.0+0356 




285.009 


3.946 


1.0 


3.6 




J1907.5+0602 


PSR b 


286.894 


6.034 


2.4 


7.4 


MGRO J1908+06 














HESS J1908+063 


J1911.0+0905 


SNR 


287.761 


9.087 


1.7 


1.5 


G43.3-0.2 


J1923.0+1411 


SNR 


290.768 


14.191 


-0.3 


3.4 


W51 














HESS J1923+141 


J1953.2+3249 


PSR 


298.325 


32.818 


-0.0 


0.0 




J1954.4+2838 


SNR 


298.614 


28.649 


0.6 


4.3 


G65.1+0.6 


J1958.1+2848 


PSR b 


299.531 


28.803 


0.1 


4.0 
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Table 1 — Continued 



Fermi LAT 








Tibet-Ill 


Milagro a 


Source 


Source 


Class 


R.A. 


Dec. 


Signi. 


Signi. 


Associations 


(OFGL) 




(dee) 


(dee) 

V 0/ 








J2001.0+4352 




300.272 


43.871 


-0.5 


-0.9 




J2020.8+3649 


PSR 


305.223 


36.830 


2.2 


12.4 


MGRO J2019+37 


J2021. 5+4026 


PSR b 


305.398 


40.439 


2.2 


4.2 




J2027.5+3334 




306.882 


33.574 


-0.3 


-0.2 




J2032. 2+4122 


PSR b 


308.058 


41.376 


2.4 


7.6 


TeV J2032+4130 














MGRO J2031+41 


J2055.5+2540 




313.895 


25.673 


-0.0 


-0.0 




J2110.8+4608 




317.702 


46.137 


0.3 


1.1 




J2214.8+3002 




333.705 


30.049 


-1.0 


0.6 




J2302. 9+4443 




345.746 


44.723 


-0.0 


-0.6 




LAT PSR J2238+59 


PSR b 


339.561 


59.080 


2.5 


4.7 





Significance of the Milagro observations. Taken from (jAbdo et al. 20 09b). 

I 

b These pulsars are newly discovered by the Fermi LAT observations (jAbdo et al. 2009a 
Abdo et al. 2009cl ). 



This puls ar is not in the Fer mi bright source list, but it is detected by the latest Fermi LAT 



observation (jAbdo et al. 2009cl ). 



